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The status of the Monte Carlo programs for the simulation of the r lepton production in high energy accelerator 
experiments and decay is reviewed. In particular, the status of the following packages is discussed: (i) TAUOLA 
for r- lepton decay, (ii) PHOTOS for radiative corrections in decays, (iii) KORALB , KORALZ , KKMC packages for r-pair 
production in e + e~ collisions and (iv) universal interface of TAUOLA for the decay of r-leptons produced by "any" 
generator. Special emphasis on requirements from new and future experiments is given. Some considerations 
about the software organization necessary to keep simultaneously distinct physics initializations for TAUOLA are 
also included. 



1. Introduction 

The TAUOLA package [Q-Q for the simulation of 
r lepton decays and PHOTOS |||] for the simu- 
lation of radiative corrections in decay, are com- 
puting projects with a rather long history. Writ- 
ten and maintained by well defined authors, they 
nonetheless migrated into a wide range of appli- 
cations where they became ingredients of com- 
plicated simulation chains. As a consequence, a 
large number of different versions are presently 
in use. From the algorithmic point of view, they 
often differ only in a few small details, but in- 
corporate a substantial amount of specific results 
from the distinct r-lepton measurements. Such 
versions were mainly maintained (and will remain 
so maintained) by the experiments taking preci- 
sion data on r leptons. On the other hand, many 
new applications were developed recently, often 
requiring program versions which differ because 
of interfaces to other packages (eg. the format of 
the event record had to be adjusted, the organi- 
zation of the input parameters changed, etc.). 

In the last two years, substantial progress for 
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the simulation of r-pair production was achieved. 
Control of the systematic errors, for the theo- 
retical predictions, at the few permille level was 
achieved (or, for some centre-fo-mass system en- 
ergies, is relatively easy to achieve in the near 
future) for the centre of mass energy from the r- 
pair production threshold up to the energies of 
the future linear colliders. 

2. Versions of TAUOLA Monte Carlo 

In ref. Q the setup for constructing specific 
versions of TAUOLA and PHOTOS from the sin- 
gle set of files was prepared and documented. 
The system was prepared for the software librar- 
ians and advanced users interested in updating 
the packages for the multipurpose environment. 
The idea was to create a repository which allows 
one to keep all main options of TAUOLA devel- 
oped for different purposes in a relatively com- 
pact form, without duplications of semi-identical 
parts. The repository was set to produce the stan- 
dard FORTRAN files which can be later handled, 
exactly the same way, as the published versions 
of the packages. 

Motivations for versioning: 

1. PHOTOS: Versions of the FORTRAN code are 
necessary because of the different versions 
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of the HEPEVT common block being in use in 
the HEP libraries (single/double precision, 
dimensionality of storage matrices). 

2. TAUOLA: Versions of the FORTRAN code are 
motivated by: (A) different versions of ini- 
tialization of physics parameters; (B) inter- 
faces to different Monte Carlo generators for 
production of r-lepton(s); and (C) different 
versions of the HEPEVT common block: 

• (A) Different physics initializations: 

(1) As published in §; 

(2) As initialized by the ALEPH col- 
laboration jTj] (it is suggested to use 
this version only with the help of the 
collaboration's advice); 

(3) As initialized by the CLEO collab- 
oration ||] (see the printout of this ver- 
sion for details); 

(4) Further coding of some individual 
decay modes. 

• (B) Different interfaces with MC gen- 
erators: 

(1) Old demo program as in the pub- 
lished version ||; 

(2) Interface to KKMC [§; 

(3) New universal interface using 
HEPEVT common block. 

• (C) Different versions of the HEPEVT 
common block. 

3. TAUOLA and PHOTOS: different versions 
of random number generators. 

4. TAUOLA and PHOTOS: makefiles with 
different compiler flags. 

Standard tools are used in the discussed setup: 
cpp the C-language pre-compiler: its if, elif 
and include commands, as well as UNIX logi- 
cal links and cat command. The aim is to pro- 
vide full backward compatibility at the level of 
the FORTRAN source with the published versions, 
or other versions being in use at present, ft 
is expected that one will use the setup to cre- 
ate her/his version of the TAUOLA and PHOTOS 
libraries and interfaces (subdirectories tauola/ 



and photos/) consisting of the plain FORTRAN 
code. The original subdirectories of the setup can 
then be erased or stored for future use. 

The following directories are created 
once the actions of the setup are completed: 

1. photos/: Standard directory with the 
FORTRAN code of the PHOTOS library and its 
demo. 

2. tauola/: Standard directory with the 
FORTRAN code of the TAUOLA library, its de- 
mos and example outputs. 

(a) tauola/demo-standalone: Demo 
program for TAUOLA executed in a 
standalone mode. 

(b) tauola/demo-jetset: Demo program 
for TAUOLA executed with universal 
interface to physics event generators 
based on the HEPEVT common block. 
In this demo HEPEVT is filled from the 
JETSET74 JkJ] Monte Carlo generator. 

(c) tauola/demo-KK-f ace: Interface to 
the KK Monte Carlo ||. 

2.1. Options for PHOTOS Monte Carlo 

The different options of PHOTOS which can be 
created^ correspond solely to the different ver- 
sions of the HEPEVT common block. The possible 
options are: 

1. KK-all - for KK Monte Carlo 

2. 2kD-all - dimension 2000 double precision 

3. 4kD-all - dimension 4000 double precision 

4. 2kR-all - dimension 2000 single precision 

5. lOkD-all - dimension 10000 double preci- 
sion 

2.2. Options for TAUOLA Monte Carlo 

Basic options for the physics initializations, 
activated with the single command make, are: 
epe; cleo; aleph. The three possible ver- 
sions of the created tauola directory correspond 

2 Executing single command: make xxx-all. 



3 



to form-factors and branching ratios defined re- 
spectively as in: (cpc) the published version of 
TAUOLA; (aleph) the conventions as adopted by 
the ALEPH collaboration, (cleo) the conventions 
as adopted by the CLEO collaboration. 

• Additional parametrizations for form-factors, 
which can be useful in some applications, are 
placed in the separate directory. At present, the 
code used in refs. and 12 1 is stored there. 

• It is often necessary to change some of 
the TAUDLA input parameters like branching ra- 
tios, the mass of the r-lepton, etc. It is 
convenient to have it done once for all appli- 
cations i.e. demo-KK-f ace, demo-jetset and 
demo-standalone. The special arrangements to 
do it, in a consistent way, are provided. For de- 
tails of the initialization routines, which are semi- 
identical in the three cases, see refs. |[|-[|. Note 
that special care from the physics point of view 
is needed. Often, input parameters are inter- 
related with the actual choice of form factors. 
The changes should be thus performed simulta- 
neously. 

3. Universal interface based on HEPEVT 
common block 

The universal interface of TAUDLA for "any" r 
production generator is provided. It uses as an 
input, the HEPEVT common block and operates 
on its content only. As a demonstration example 
the interface is combined with the JETSET gener- 
ator, however it should work in the same manner 
with the PYTHIA0, HERWIG or ISAJET generators 
as well. 

The interface work in the following way: 

• The T-lepton should be forced to be stable 
in the r production generator. 

• The content of the HEPEVT common block 
is searched for all r leptons and neutrinos 
first. 

• It is checked if there are r-flavour pairs (two 
r-leptons or r-lepton and r-neutrino) orig- 
inating from the same mother. 



• The decays of the r-flavour pairs are per- 
formed with TAUOLA. Longitudinal spin ef- 
fects are generated in the case of the r pro- 
duced from decay of: W — > tv, Z/7 — > rr, 
the neutral Higgs boson H — > rr, and the 



charged Higgs boson H d 



tv. Parallel 
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or anti-parallel spin configurations are gen- 
erated, before calling on the r decay, and 
then the decays of 100 % polarized r's are 
executed. 

• In the case of the Higgs boson (for the spin 
correlations to be generated) the identifier 
of the r mother must be that of the Higgs. 

• In case of the W and Z/7 it is not neces- 
sary. If from the same mother as that of 
the r there is produced also a v T , the W is 
reconstructed by the interface as the sum of 
the two. Simmilarily, Z/7 is reconstructed 
if from the same mother another r is pro- 
duced. 

• Let us note that the calculation of the r po- 
larization created from the Z and/or virtual 
7 (as function of the direction) represents 
a rather non-trivial extension. Dedicated 
study of the production matrix elements of 
the host generator is necessary. A separate 
paper |l| is devoted to this point; let us 
show however some preliminary results in 
the next section. At this moment the dis- 
tribution version of the package, available 
from the www-page, does not include this 
possibility. The updated version can be ob- 
tained from the authors upon individual re- 
quests only. 

• Photon radiation in decay is performed with 
PHOTOS. 

4. Tau polarization from the Z/7 mecha- 
nism 

The best way of calculating spin state of any 
final state is with the help of the matrix element 
and the rigorous density matrix treatment. This 
is however not always possible or necessary. Of- 
ten, like in the case of the production and decay 
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of particles in the ultra-relativistic limit a simpli- 
fied approach can be sufficient. Such an approach 
was developed for KORALZ Monte Carlo program 
Hj| and its limitations were studied with the help 
of matrix element calculations of order a Jlgl . In 
fl3|l we study the question of whether a similar 
approach can not be generalized, and the approx- 
imate spin correlation calculated from the infor- 
mation stored in the HEPEVT common block by 
the typical r production program: 

The approximation consists of calculat- 
ing/approximating variables and later using the 
spin correlation of the elementary 2^2 body 
process e + e~(qq) — > t + t~, buried inside multi- 
body production process. Let us stress that such 
a procedure can never be fully controlled as its 
uncertainties and even functioning depends on 
the way the particular production program fills 
the HEPEVT common block. It will be always 
responsibility of the user to check if in the par- 
ticular case the approximation can be useful. 
Nonetheless our aim is not to replace the ma- 
trix element calculations, but rather to provide 
a method of calculating/estimating spin effects 
in cases when otherwise spin effects would not 
be taken care of, at all. Needless to say such an 
approach is limited for the spin treatment to the 
approximation not better than the leading-log. 

The principle of calculating kinematical vari- 
ables is simple. The 4-momenta of the 2-^2 
body process have to be found. The 4-momenta of 
the outcoming r's are used directly. Initial state 
momenta are constructed from the incoming and 
outcoming momenta of the particles (or fields) 
accompanying production of the Z/j state^. We 
group them accordingly to fermion number flow, 
and ambiguous additional particles are grouped 
(summed) into effective quarks to minimize their 
virtualities. Such an approach is consistent in the 
case of emission of photons or gluons with the 
leading log approximation. 

In the following, let us give a few examples of 
tests we have performed so far on the elements of 
the calculation, see ref |13l for more details. In 



fig. 1 we plot the t polarization in 2 — > 2 body 
scattering as a function of the scattering angle 
respectively for e~e + , UU, dd initial states. In 
fig. 2 we use the most spin sensitive decay mode 
t — > istt and we show the ir energy spectrum (in 
the Z rest- frame), and energy-energy correlations 
between ir + and ir~ (also in the Z rest- frame). 
Due to the Z interaction the r leptons are polar- 
ized; there is thus a shift toward lower energies 
in the n energy spectrum. From the second part 
of the plot, we see that the configurations where 
both 7r + and 7r~ have energies bigger or smaller 
than average are more populous than the mixed 
configurations, exhibiting the vector nature of the 
Z/j — t — t vertex. On the contrary, in the case of 
the production from the Higgs meson, the oppo- 
site would be true, and the mixed fast-slow (slow- 
fast) configurations would be more common, the 
energy spectrum of the it would be flat. 

5. The t pair production in e + e~ collisions 

Until 1999 the two widely used programs for 
the r-pair production in e + e~ collisions were 
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KQRALB (16| and KORALZ |lj . The first of the two 
programs was used mainly at lower energies where 
its full treatment of spin and T-mass effects was 
more important than its limited, first-order only, 
treatment of radiative corrections. KORALZ was 
more appropriate for the higher energy zone, es- 
pecially on the peak of the Z resonance, where the 
ultra-relativistic approximation of the spin effects 
was sufficient, effects of the higher order QED 
corrections were important, and, thanks to the Z 
life-time, the interference of initial and final state 
bremsstrahlung was strongly suppressed. This 
second assumption turned out not to withstand 
the numerical requirements of the LEP2 physics 
see eg. 0. Also, it is not enough for require- 
ments of the future high energy Linear Colliders 
requirements such as TESLA JL8| . 

In the recent years, a new Monte Carlo pro- 
gram, based on exponentiation performed at the 
spin amplitude level, was developed g]. Thanks 
to the applied technique, it can provide predic- 
tions at the precision level of a few permille, for 
centre-of-mass energies from r production thresh- 
old up to the energies of the linear collider range. 
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Pi gure 1. Tests of the TAU0LA universal interface. The 
r lepton polarization as a function of cos#. We have 
used y/s = Mz, M z ± 3 GeV. 



So far, precision tests were performed for the 
centre-of-mass energies corresponding to LEP2. 
For other energies, points such as matching vir- 
tual pair corrections with simulations of 4-fcrmion 



final states jlfj , photonic vacuum polarization ef- 
fects at low energies, Coulomb interactions, mass 
terms in some virtual corrections, etc. have to be 
discussed before precision at the permille level can 
be granted. For more details and comparisons be- 
tween KKMC, KORALZ and K0RALB, see table 1 and 
ref. fl. 
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Feature 


KORALB 


KORALZ 


KK 4.13 


KK 2000+? 


QED type 


O(a) 


EEX 


CEEX, EEX 


CEEX, EEX 


CEEX(ISR+FSR) 


none 


none 


{a, ah; c* 2 L 2 , o 2 !, 1 } 


{...a 2 il it> 3 I, 3 } 


EEX(ISR'FSR) 


none 




{ C ,,aL,a 2 L 2 , C> 3 L 3 } 


{...o 2 L 2 ,a 3 L 3 } 


ISR-FSR int. 


0(a) 




qL > C EEX 


{a, aL } CEEX 




1 7 


1 , 2coll. 7 


1,2, 3coll. 7 


up to 3 7 


Electroweak 


No 2-res. 


DIZET 6.x 


DIZET 6.x 


New version? 


Beam polar. 


long+trans. 


longit. 


long+trans. 


long+trans. 


r polar. 




longit. 


long+trans. 


long+trans. 


Hadronization 




JETSET 


JETSET 


PYTHIA 


t decay 


TAUOLA 


TAUOLA 


TAUOLA 


TAUOLA 






No 


Yes 


Yes 


Beamstrahlung 




No 


Yes 


Yes 


Beam spread 




No 


Yes 


Yes 






Yes 


No 


Yes 


ee channel 




No 


No 


Yes 


f f channel 




No 


No 


yes? 


WW channel 




No 


No 


yes? 



Table 1: Comparison of the main features for the 
KORALB, KORALZ, and KKMC programs, see ref. [§ for 
more details. 



6. Summary and future possibilities 

The status of the programs for the production 
of t leptons was reviewed. The high precision 
Monte Carlo program KKMC for r-pair production 
in e + e~ collisions was presented. The general 
purpose interface, with the partial treatment of 
spin effects, of TAUOLA to any production program 
was also reviewed. At present, in its full content, 
it is available from the authors upon individual 
request only. 

Distinct versions of the TAUOLA library for r 
lepton decay, and of PHOTOS for radiative correc- 
tions in decays, are now in use. The versions 
differ either due to physics initialization or re- 
quirements of the interfaces. We have presented 
the system for creating the required version of 
the PHOTOS and TAUOLA packages from the sin- 
gle master copies. The master copies are kept in 
relatively compact and clear form, without code 
duplications and with the help of the epp pre- 
compiler. 

We expect that this organization, on top of 
the practical goals of the every day applications, 
may serve also as a step to develop the packages 
(without loss of their physics content) into new in- 
ternal architectures. Some additional experience 



was collected in ]20|] in context of 00 C++. We 
have found that the language translation for the 
given program version is relatively easy. On the 
contrary, the question of project continuity into 
further upgrades motivated by the physics needs 
to be thought over carefully. Matching the pro- 
gramming styles, e.g. of the 00 C++ experts, with 
the strategies of testing the numerical correctness 
of consecutive versions is a rather crucial issue 
which has to be addressed. Tools and methods 
embodied in FORTRAN survive such translations 
with difficulty. 

The successful strategy will probably require 
simultaneous fluency, at a certain step of the 
project development, in FORTRAN, 00 languages 
and the physics content of the project, by all in- 
volved persons. Platform independent tools for 
mixing code in FORTRAN and 00 languages will be 
of great help. 
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